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Abstract. In the context of f(R) theories of gravity, we study the cosmological evolution of 
scalar perturbations by using a completely general procedure. We find that the exact fourth-order 
differential equation for the matter density perturbations in the longitudinal gauge, reduces to a 
second-order equation for sub-Hubble modes. This simplification is compared with the standard 
(quasi-static) equation used in the literature. We show that for general f(R) functions the quasi- 
static approximation is not justified. However for those f(R) adequately describing the present 
phase of accelerated expansion and satisfying local gravity tests, it does give a correct description 
for the evolution of perturbations. 
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INTRODUCTION 

The present phase of accelerated expansion of the universe [1] poses one of the most 
important problems of modern cosmology. It is well known that ordinary Einstein's 
equations in either a matter or radiation dominated universe give rise to decelerated 
periods of expansion. In order to have acceleration, the total energy-momentum tensor 
appearing on the right hand side of the equations should be dominated at late times by a 
hypothetical negative pressure fluid usually called dark energy fjfa. 

However, there are other possibilities to generate a period of acceleration where 
Einstein's gravity itself is modified. In one of such possibilities, new functions of the 
scalar curvature (f(R) terms) are included in the gravitational action. Although such 
theories are able to describe the accelerated expansion on cosmological scales correctly, 
they typically give rise to strong effects on smaller scales, but nevertheless, viable 
models can be constructed imposing strong constraints over f(R) [H]: /rr > for high 
curvatures, 1 + /r > for all R, /r < to ensure that ordinary General Relativity 
behaviour is recovered at early times and \/r\ <C 1 at recent epochs to ensure that local 
gravity tests hold, /r and /rr mean first and second derivative respectively of f(R) with 
respect to the argument R. 

The important question that arises is therefore how to discriminate dark energy mod- 
els from modified gravities using observations. It is known that by choosing particular 
f(R) functions, one can mimic any expansion history, and in particular that of ACDM. 
Accordingly, the exclusive use of observations from SNIa fl], baryon acoustic oscilla- 
tions [0] or CMB shift factor [5], based on different distance measurements which are 
sensitive only to the expansion history, cannot settle the question of the nature of dark 
energy [0]. 



However, there exists a different type of observations which are sensitive, not only to 
the expansion history, but also to the evolution of matter density perturbations. The fact 
that the evolution of perturbations depends on the specific gravity model, i.e. it differs 
in general from that of Einstein's gravity even though the background evolution is the 
same, means that this kind of observations will help distinguishing between different 
models for acceleration. 

In this work [|7Q we find the exact equation for the evolution of matter density pertur- 
bations for arbitrary f(R) theories. Such problem had been previously considered in the 
literature [8] and approximated equations have been widely used based on the so called 
quasi- static approximation in which all the time derivative terms for the gravitational po- 
tentials are discarded, and only those including density perturbations are kept. From our 
exact result, we will be able to determine under which conditions such an approximation 
can be justified. 



DENSITY PERTURBATIONS EVOLUTION IN f(R) THEORIES 



The corresponding equations of motion for that action are modified with respect to their 
usual General Relativity counterpartners. In our research we use a flat Robertson- Walker 
perturbed metric in the longitudinal gauge 



where 4> = <J>(t] , it) and \F = *P(tj , it) are the scalar perturbations. Using this perturbed 
metric and a perturbed energy-momentum tensor for dust matter perfect fluid, the first 
order perturbed equations may be written for a general f(R) function assuming that 
background equations hold. A relevant point of these theories is that potentials <J» and 
*F are not equal provided /rr ^ 0. From Einstein's equations, it is possible to derive a 
fourth order differential equation for matter density perturbation 8 = 8p/po - with po 
unperturbed density and 8p density constrast with respect to the background -alone. The 
process can be found in y\\. The resulting equation can be written as follows: 



J 84 )/ 5 iv + J 8 3)/ 5 , "+(a 2 + j8 2 , / )5" + (a 1 + J 8 1)/ )5'+(ao + j8o ) /)5 = (3) 



where the coefficients jS^/ (i = 1, ...,4) involve terms with f R and f' R . Equivalently, 0f z 
(z = 0, 1,2) contain terms coming from the linear part (Einstein-Hilbert) of f(R) in R 
(background curvature). Parameter e = Jif/k = a' / (ak), where ' denotes derivative with 
respect to conformal time 7], will be useful to perform a perturbative expansion of the 
previous coefficients oc's and /3's in the sub-Hubble limit. We also define the following 



dimensionless parameters: Kj = / Jf 1+1 and fi = f' R {J) / (Jtfjfit) (i, j = 1,2, 3). 



Let us consider the modified gravitational action: 




(1) 



ds 2 = a 2 (ri ) [(1+ 24>) Jt] 2 -(1- 2¥) {dr 2 + r 2 dQ.l)] 



(2) 



SUB-HUBBLE MODES AND THE QUASI-STATIC LIMIT 



We are interested in the possible effects on the growth of density perturbations once they 
enter the Hubble radius in the matter dominated era. In the sub-Hubble limit £ <C 1, it 
can be seen that the fa j and fo f coefficients are supressed by £ 2 with respect to /32,/, 
J3i,/ and /3o.f, i.e., in this limit the equation for perturbations reduces to the following 
second order expression: 



s ;z , WK< i (l+fR) 5 ^ 2 (-l + ^)(^-K2)-$f 4 RR (K2-2)k*8nGp a\ 

(l+/,) 5 (-l + K 1 ) + ^(l+^)(^2-2)^ d " U(4) 

where we have taken only the leading terms in the £ expansion for the a and /3 
coefficients. This expression can be compared with that usually considered in literature, 
obtained after performing strong simplifications in the perturbed equations by neglecting 
time derivatives of <J» and *F potentials. Thus in D9D they obtain: 

a L I+JR 

We will now restrict ourselves to models satisfying all the viability conditions mentioned 
in the introduction, including |/^| <C 1. 

Note that the quasi-static expression © is only recovered in the matter era (i.e. for 
Jt? = 2/rj) or for a pure ACDM evolution for the background dynamics. Nevertheless 
in the limit | fa |<C 1 it can be proven using the background equations of motion that 
1 + Ki — K2 ~ and therefore 2k\ — ~ —2 + K2 ~ —1 + K\ what is nothing but the 
fact that for viable models the background evolution resembles that of ACDM. This this 
fact allows to simplify expression © to approximately become ©. 

In other words, although for general f(R) functions the quasi-static approximation is 
not justified, for those viable functions describing the present phase of accelerated ex- 
pansion and satisfying local gravity tests, it gives a correct description for the evolution 
of perturbations. 

In order to check our results we proposed two particular f(R) models to obtain solu- 
tions both using © and ©. As commented before, for viable models the background 
evolution resemble that of ACDM at low redshifts and that of a matter dominated uni- 
verse at high redshifts. Nevertheless the f(R) contribution gives the dominant contribu- 
tion for small curvatures and therefore it may explain the cosmological acceleration. 

The two studied models are of the form f(R) == c\R p (where units Hq will be 
considered). According to the results presented in [10], models of this type include both 
matter dominated and late-time accelerated universe provided the parameters satisfy 
c\ <0and0<p< 1. 

The first model (left figure) will have parameters c\ = —4 and p = 0.63 and it is not 
accomplishing condition | fn |<C 1. Second model (right figure) will have c\ = —4.3 
and p = 0.01 and it does accomplish condition | /r |<C 1. For each model, we have 
compared our result with the standard ACDM and the quasi-static approximation Q. 
We saw that for first model quasi-static and correct approximation do not give the 



same evolution but for the second one the quasi- static approximation gives a correct 
description for the evolution. 
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Left figure: 8^ with k — 0.2hMpc~ l for first model and ACDM. Both, standard quasi-static evolution 
and equation © have been plotted and they clearly differ. ; Right figure: 8^ with k = 1. 67 hMpc~ l for 
second model and ACDM. The quasi-static evolution is indistinguishable from that coming from ©, but 
both diverge from ACDM behaviour as z decreases. 

CONCLUSIONS 

We have shown that for sub-Hubble modes, the differential equation for the evolution of 
density perturbation reduces to a second order equation and compared this result with 
that obtained within quasi-static approximation used in the literature and found that for 
arbitrary f(R) functions, such an approximation is not justified. 

However for theories with <C 1 today, perturbative calculation for sub-Hubble 
modes requires to take into account, not only the first terms, but also higher-order terms 
in e = J$?/k. In that case, the resummation of such terms modifies the equation which 
can be seen to be equivalent to the quasi-static case, but only if the universe expands as 
in a matter dominated phase or in a ACDM model. Finally, the fact that for models with 
\f R \ <C 1 the background behaves today precisely as that of ACDM makes the quasi- 
static approximation correct in those cases. 
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